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O Abstract 

The purpose of this paper is to propose an inexpensive micromechanics-based scheme 
for stiffness homogenization and strength estimation of mortars containing crushed bricks, 
known as cocciopesto. The model utiHzes the Mori-Tanaka method to determine the effective 

I— I stiffness, combined with estimates of quadratic invariants of the deviatoric stresses inside 

* ^ phases to predict the compressive strength. Special attention is paid to the representation 

of C-S-H gel layer around bricks and interfacial transition zone around sand aggregates, 

1— H which renders the predictions sensitive to particle sizes. Several parametric studies are 

performed to demonstrate that the method correctly reproduces data and trends reported 
^ in available literature. Moreover, the model is based exclusively on parameters with clear 

^ physical or geometrical meaning and as such it provides a convenient framework for its further 

experimental validation. 

e 

Keywords, micromechanics, homogenization, strength estimation, cocciopesto, C-S-H gel 
C coating, iterfacial transition zone 

o 

1 Introduction 

(N 

The use of lime as a binder in mortars is associated with well-known inconveniences such as slow 
l/^ setting and carbonation, high drying shrinkage and porosity and low mechanical strength [1]. 

Although these limitations have been overcome with the use of Portland cement in the last 
50 years, Hme mortars still find use in the restoration of historic structures. This is mainly due 
to their superior compatibility with the original materials in contrast to many modern renovation 
^ render systems, e.g. [14, 15, 27]. 

Mechanical properties of lime mortars can be improved by a suitable design of the mixture. 
>■ Phoenicians were probably the first ones who added crushed clay products, such as burnt bricks, 

tiles or pieces of pottery, to lime mortars to increase their durability and strength. Romans 
called such material cocciopesto and utilized this mortar in areas where other natural pozzolans 



were not available. The cocciopesto-based structures exhibit increased ductility, leading to their 
remarkable resistance to earthquakes [3, 19]. 

Much later, it was found that the mortars containing crushed clay bricks, burnt at 600-900°C, 
exhibit a hydraulic character, manifested by the formation of a thin layer of Calcium-Silicate- 
Hydrate (C-S-H) gel at the lime-brick interface [18]. Since C-S-H gel is the key component 
responsible for favorable mechanical performance of Portland cement pastes [21], it is generally 
conjectured that the enhanced performance of cocciopesto mortars can be attributed to the high 
strength and stiffness of the C-S-H gel coating [3, 18, 19, 27]. This mechanism competes with 
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Figure 1: Scheme of the micromechanics-based model. The numbers in parentheses refer to 
indexes of individual phases. 

the formation of the Interfacial Transition Zone (ITZ) at the matrix-aggregate interface, known 
to possess higher porosity and thus lower stiffness in cement-based mortars, e.g. [22, 26, 34]. 

The purpose of this work is to interpret these experimental findings by a micromechan- 
ical model based on the Mori-Tanaka method [17], motivated by its recent applications to 
related material systems. These include, for example, estimates of effective thermal conduc- 
tivity of rubber-reinforced cement composites [29] , elasticity predictions for early-age cement [5] 
or alkali-activated [33] pastes, upscaling the compressive strength of cement mortars [24], and 
multi-scale simulations of three-point bending tests of concrete specimens [32]. Here, we exploit 
these developments to propose a simple analytical model for stiffness and strength estimation 
of cocciopesto mortars in Section 2. In particular, the elasticity predictions utilize Benveniste's 
reformulation [4] of the Mori-Tanaka method [17], whereas the strength predictions build on re- 
cent results by Pichler and Hellmich [24] , who demonstrated that compressive strength is closely 
related to the quadratic average of the deviatoric stress in the weakest phase. A particular at- 
tention is paid to representation of coatings by C-S-H gel and ITZ, which renders the predictions 
sensitive to the size of brick particles and aggregates. In Section 3, we verify predictions of the 
proposed scheme against data available in open literature. These finding are summarized in 
Section 4, mainly as a support for future validation of the model against experimental results. 
Finally, in Appendix A we gather technical details needed to account for coated inclusions in 
order to make the paper self-contained. 

In what follows, the Mandel representation of symmetric tensorial quantities is systematically 
employed, e.g. [16, p. 23]. In particular, italic letters, e.g. a or A, refer to scalar quantities and 
boldface letters, e.g. a or A, denote vectors or matrix representations of second- or fourth-order 
tensors. and (A)^^ standardly denote the matrix transpose and the inverse matrix. Other 
symbols and abbreviations are introduced in the text when needed. 



We consider a composite sample occupying domain ri, composed of n distinct phases indexed 
by r. The value r = is reserved for the matrix phase and r = 1, . . . , n refer to heterogeneities 
having the shape of a sphere or spherical shell, see Fig. 1. Volume fraction of the r-th phase are 



2 Model 




Elastic properties of individual phases are specified by the material stiffness matrix L^'^\ As 
each phase is assumed to be homogeneous and isotropic, we have 




where K^'^^ and G^^^ are the bulk and shear moduli of the r-th phase, and ly and Id denote 
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the orthogonal projections to the volumetric and deviatoric components, e.g. [12], so that 



e{x) = {ly + Ib)£{x) = ey{x)l + ei){x), (2a) 
ct{x) = {ly + Iy))(t{x) = fTv(a;)l + ctj:,{x), (2b) 

for X ^ n. In Eq. (2), e and cr refer to local stresses and strains, £y and £d are the volumetric 
and deviatoric strain components, ay and ctd refer to the stress components and 1 is the second- 
order unit tensor (in the matrix representation). 

Development of the model follows the standard routine of the continuum micromechanics, 
e.g. [35]. The sample O is subjected to the overall strain loading E. Neglecting the interaction 
among phases, the mean strains inside heterogeneities are obtained as 

i;M = Ai^}£; forr = l,...,n, 

(r) 

where A^jj is the dilute concentration factor of the r-th phase, see Section 2.1. In Section 2.2, 
after accounting for the phase interaction, these are combined to the full concentration factors 
satisfying 

^(r)^^(r)^ for r = 0,...,n, 

utilized next to estimate the overall stiffness of the composite material, L^s- Moreover, as 
outlined in Section 2.3, expression for the overall stiffness also encodes the mean value of the 
quadratic invariant of the local stress deviator ctd defined as 



4^ = ^^J^^^^cr^{xya^{x)dx, (3) 

that can be directly used to estimate the overall strength of a material. 
2.1 Dilute concentration factors 

Due to geometrical and material isotropy of individual phases, the dilute concentration factors 
attain the form analogous to (1): 

^diJ = ^Sv^v + ^Sd^d for r = 1, . . . , n. (4) 

The expressions for the components are given separately for the uncoated (r = 1) and coated (r = 
2, . . . , 5) particles. Namely, in the first case it holds 



A 



(1) 

dil,V 



^dii,D Q(o) ^p{o)^Q(i) _Q(o)y 

with the auxiliary factors following from the Eshelby solution [9] in the form 

«(0) = 1 + ^-^°^ «(0) ^ 2(4-5.W) 

3(1 + 1/(0))' P 15(1-1/(0))' 

where zy(0) is the Poisson ratio of the matrix phase. 

The coated case is more involved, and was first solved in its full generality by Herve and 
Zaoui [10] for multi-layered spherical inclusion. To apply their results in the current setting, we 
locally number the phases by the index i = [ii, i2, is]"*", see Fig. 2, where i = [2,3,0]"'' for the 
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Figure 2: Scheme of a single-layer inclusion 



Table 1: Reference properties of individual phases; p denotes density, ft is tensile 
strength, m is the mass fraction. 



r 


Phase 


P 


E V ft 


m 


R Note 






[kgm^^] 


[MPa] [-] [MPa] 


[-] 


[\im\ 





Pure lime matrix 


1,200^ 


2, 000^^1 0.25^81 0.4I8J 


3 


X 


1 


Voids 


X 


10-9 0.25 X 


X 


[20] ^(1) ^ 35%a 


2 


Clay brick 


2, 300^ 


5,000^ 0.17^ 3.2^ 


1 


500 


3 


C-S-H gel 


2,000[ii] 22,000[^] 0.2['^1 x 


X 


510[61 


4 


Siliceous sand 


2, 700^ 


60,000[32l0.17[32l 48P11 


1 


500 


5 


ITZ 


1,200'' 


500= 0.25'' X 


X 


520^= 



^ Our own (unpublished) data. 

^ Same value as for the lime matrix. 

Set as in [34] for cement-based concretes, i.e. Young's modulus to 20-40% of the 

value for the matrix phase and thickness to 20 [xm. 



brick-C-S-H conglomerate and i = [4, 5, 0]""" refers to a sand particle coated by ITZ. Now, we 
have 



^dil,V ~ n2 ' ^dil,V ~ r)2 ' ^'-'l 

^11 Vll 



)1 



and 

i(n) 



^^ = A,--^-^B,. (6) 

r \ 21 R(«2)5 _ E>(i2-1)5 

^dil,D ^2 c . n,,('ion^D('io~),'? D^io-lLSA 2, 



where the auxiliary factors are provided in Appendix A. 



2.2 Stiffness estimates 

In Benveniste's [4] interpretation of the original Mori-Tanaka method [17], the mutual interaction 
among heterogeneities is modeled by loading each particle by the average strain in the matrix 
phase E^^^ instead of E. For this purpose, we relate E^^^ to -E by a strain compatibility 
condition, valid under the dilute approximation, 

n 
r=l 
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from which we express the full concentration factors as 

n 

A(0) = (e(0)/ + ^eMAS)-\ 

r=l 

= aJ:;a(0) for r = l,...,n. 

Utilizing a universal relation 

n 

Leff = 5^cMi:WAW, 

r=0 

we can see that the effective stiffness inherits the symmetry of individual phases (1) with 

^eff = r\ ) 

c(0) + E".icM4;,v 

Gcff = —■ (7b) 

2.3 Strength estimates 

As recognized first by Kreher [13], the fluctuations of stresses and strains in individual phases 
can be estimated from the energy conservation condition due 



+ / 2G^'-hl{x)eiy{x)dx), 



(8) 



expressing the conservation of energy on macroscale due to E and the average local values due 
to ev s-iicl £d- Differentiating (8) with respect to G^^^ ^ we obtain 



for r = 0, . . . , n. Next, we recognize that a^j^^ = {l/2G^'''^)e''^^ and recall the definition of the 
quadratic invariant (3) to arrive at 



^ ^ 



As thoroughly demonstrated by and Pichler and Hellmich [24], this quantity is closely 
related to the compressive strength /c of cement pastes at various degrees of hydration. Here, 
we postulate that 

UPi) Jt\p2) .... 

where w = 0, . . . ,n is the index of the weakest phase and p refers to a parameter characterizing 
the mixture composition, see the next section for concrete examples. 
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3 Results and discussion 



The purpose of this section is to examine the trends in mechanical properties as predicted by the 
proposed scheme. Default data for individual phases, summarized in Table 1, were partly as- 
sembled from open literature and complemented with our own, yet unpublished, measurements. 
Note that the matrix-brick-sand fractions correspond to a typical composition of historic lime 
mortars [3, 2] and that the engineering constants E and v are connected to the bulk and shear 
moduli through well-known relations, e.g. [16, p. 23], 



3{l-2u)' 2(1 + 1/)' 

, motivated by experimental findings in e.g. [3, 31, 30], we assume that an increase in the C- 
S-H gel volume fraction (Ac^^^) is compensated by the corresponding changes for matrix (Ac^^^), 
voids (Ac^^)), and clay bricks (Ac*^^^), so that 

Ac(o) + Ac(i) + Ac(2) + Ac(3) = 0, 

where we set for simplicity Ac^^'^ = Ac^^^ = Ac^^\ Analogously, an increase in the volume 
fraction of ITZ corresponds to the decrease in volume of the matrix phase: 

Ac(o) + Ac(5) = 0. 

In the strength estimates, the imposed loading simulates the uniaxial compression test, for 
which E = [-1, 0, 0, 0, 0, 0]"^ and the average strain follows from 

E = (Leff)-^S. 

We assume that ITZ is the weakest phase, i.e. w = 5 m Eq. (10), and estimate the derivative in 
Eq. (9) by the forward difference with the step size of AG(5) = 1 Fa} 

3.1 Effect of coatings 

The first aspect we would like to discuss is the effect of coating on the dilute concentration factors 
of the brick and sand particles. Fig. 3 demonstrates that, in terms of the volumetric phase strains, 
the effects of C-S-H and ITZ are comparable, despite the fact the C-S-H is stiffer and ITZ more 
complaint than the matrix phase. The differences in the deviatoric part, which drive the strength 
estimates according to Eq. (10), become more pronounced with the increasing thickness. This 
indicates that the contribution of brick and sand particles to the overall properties might still 
be different, once accounting phase properties and their interaction, see Section 3.3 for further 
discussion. 

3.2 Influence of porosity 

Analogously to the cement pastes, porosity has a major influence on the overall properties of 
lime-based mortars. This is also confirmed by the results of the proposed model shown in Fig. 4. 
As for the overall stiffness, for the realistic range of porosities of 25-40% [8], the estimates (7) 
predict the values of Young modulus between ~ 2,000 and 1,000 MPa. This is consistent with 
the values reported in [3] historic lime mortars (without pozzolan admixtures). 

As for the strength estimates, it follows from Fig. 4 (b) that they reproduce [21, p. 280] 

/c(c«)=/c(0)(l-cW)", (11) 

^ Our results are reproducible with a MATLAB code Homogenizator MT, freely available at http://mech.fsv. 
cvut . cz/~nezerka/ software. 
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Influence of the coating thickness on the dilute concentration factors for (a) brick and 
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Figure 4: (a) Stiffness- and (b) strength-porosity relations. 



with n = 1.04, yielding practically the linear strength-porosity scaling. Unfortunately, we are 
currently unable to validate this prediction against experiments; the only available work we 
are aware of by Papayianni and Stefanidou [23] does not contain enough data. Still, Eq. (11) 
complies with the fact the influence of porosity is much smaller in lime mortars than in cement- 
based materials , for which 3 is typically used, . 



3.3 Size effects 

Now we proceed to clarify the impact of brick and sand particles on the overall mechanical 
properties.^ In particular, when increasing the size of brick particles, the material becomes 
more complaint since the stiffening effect of C-S-H layer decreases. Fig. 5(a). This also increases 
the deviatoric stresses in ITZ, as manifested by the strength reduction visible in Fig. 5(b). 
These effects practically stabilize for particles larger than 0.5 mm and their magnitude is rather 
limited: the stiffness decreases by about 10 % and the strength only by 4 %. Such trends are 
qualitatively consistent with the results presented e.g. in [3, 19, 27]. 

Larger sand particles, on the other hand, tend to make the composite material stiffer. 
Fig. 6(a), by compensating for inferior mechanical properties of ITZ. Since the relative thickness 
of ITZ layer decreases, the stresses inside this phase increase and the material becomes weaker in 
overall. Fig. 6(b). When compared to brick particles, these effects are much more pronounced: 
in the considered range of radii, the Young modulus increases by about 100 % and the strength 
decreases by 25 % with no tendency to stabilize. This agrees well with experimental outcomes 
reported in [28]. 



4 Conclusions 

In the present work, following the recent developments presented in [25, 24, 33, 32], a simple 
micromechanics-based scheme for strength and stiffness estimates of cocciopesto mortars has 
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been presented. The model utilizes directly measurable material and geometrical properties of 
individual phases and is free of adjustable parameters. On the basis of presented results, we 
conclude that the model 

1. predicts realistic values of the overall Young modulus and the strength-porosity scaling, 

2. captures the "smaller is stiffer" and "smaller is stronger" trends for crushed brick particles, 

3. captures the "larger is stiffer" and "larger is weaker" trends for sand aggregates, 

4. explains the positive role of crushed bricks in comparison with sand aggregates. 

Of course, in order to accept this model for practical use, it needs to be validated against 
comprehensive experimental data at micro- and macro-scales. This topic is currently under 
investigation and will be reported separately. 



Acknowledgments This work was supported by the Grant Agency of the Czech Technical 
University in Prague, project No. SGS12/027/OHK1/1T/11 (VN, JZ) and by the Ministry of 
Education, Youth and Sports of the Czech Republic, project No. 684077003 (JZ). 



References 

[1] A. Arizzi and G. Cultrone. Aerial lime-based mortars blended with a pozzolanic additive 
and different admixtures: A miner alogical, textural and physical-mechanical study. Con- 
struction and Building Materials, 31:135-143, 2012. 

[2] G. Baronio and L. Binda. Study of the pozzolanicity of some bricks and clays. Construction 
and Building Materials, 11:41-46, 1997. 

[3] G. Baronio, L. Binda, and N. Lombardini. The role of brick pebbles and dust in conglom- 
erates based on hydrated lime and crushed bricks. Construction and Building Materials, 
11:33-40, 1997. 



8 



[4] Y. Benveniste. A new approach to the apphcation of Mori-Tanaka theory in composite 
materials. Mechanics of Materials, 6:147-157, 1987. 

[5] O. Bernard, Ulm F.-J., and E. Lemarchand. A multiscale micromechanics-hydration model 
for the early-age elastic properties of cement-based materials. Cement and Concrete Re- 
search, 33:1293-1309, 2003. 

[6] H. Bokc, S. Akkurt, B. Ipekoglu, and E. Ugurlu. Characteristics of brick used as aggregate 
in historic brick-lime mortars and plasters. Cement and Concrete Research, 36:1115-1122, 
2006. 

[7] G. Constantinides and F.J. Ulm. The effect of two types of C-S-H on the elasticity of cement- 
based materials: Results from nanoindentation and micromechanical modeling. Cement and 
Concrete Research, 34:67-80, 2004. 

[8] M. F. Drdacky and D. Michoinova. Lime mortars with natural fibres. In A. M. Brandt, 

V. C. Li, and L H. Marshall, editors, Brittle Matrix Composites 7: Proceedings of the 7th 
International Symposium, pages 523-532, Cambridge, UK, 2003. Institute of Fundamental 
Technological Research, Woodhead Publishing Limited. 

[9] J.D. Eshelby. The Determination of the Elastic Field of an Ellipsoidal Inclusion, and Related 
Problems. Proceedings of the Royal Society of London, 241:376-396, 1957. 

[10] E. Herve and A. Zaoui. n-layered inclusion-based micromechanical modelling. International 
Journal of Engineering Science, 31:1-10, 1993. 

[11] J.T. Jeffrey and M.J. Hamlin. A colloidal interpretation of chemical aging of the C-S-H gel 
and its effects on the properties of cement paste. Cement and Concrete Research, 36:30-38, 
2006. 

[12] M. Jirasek. Basic concepts and equations of solid mechanics. Revue Europeenne de Genie 
Civil, 11:879-892, 2007. 

[13] W. Kreher. Residual stresses and stored elastic energy of composites and polycrystals. 
Journal of the Mechanics and Physics of Solids, 38:115-128, 1990. 

[14] P. Maravelaki-Kalaitzaki, A. Bakolas, and A. Moropoulou. Physico-chemical study of Cre- 
tan ancient mortars. Cement and Concrete Research, 33:651-661, 2003. 

[15] D. Michoniova. Questions about renovation plasters. Zprdvy pamdtkove pece, 65:313-316, 
2005. in Czech. 

[16] G. W. Milton. The Theory of Composites. Cambridge Monographs on Applied and Com- 
putational Mathematics. Cambridge University Press, 2002. 

[17] T. Mori and Tanaka K. Average stress in matrix and average elastic energy of materials 
with mixfitting inclusions. Acta Metallurgica, 21:571-574, 1973. 

[18] A. Moropoulou, A. Bakolas, and K. Bisbikou. Characterization of ancient, Byzantine and 
later historic mortars by thermal and X-ray diffraction techniques. Thermochemica Acta, 
269/270:779-995, 1995. 

[19] A. Moropoulou, A.S. Cakmak, G. Biscontin, A. Bakolas, and E. Zendri. Advanced Byzan- 
tine cement based composites resisting earthquake stresses: the crushed brick/lime mortars 
of Justinian's Hagia Sophia. Construction and Building Materials, 16(8) :543 - 552, 2002. 



9 



[20] M.J. Mosquera, B. Silva, B. Pricto, and E. Ruiz-Herrera. Addition of cement to lime-based 
mortars: Effect on pore structure and vapor transport. Cement and Concrete Research, 
36:1635-1642, 2006. 

[21] A.M. Neville. Properties of Concrete: Fourth Edition. Longman, 1996. 

[22] J. P. Ollivier, J. C. Maso, and B. Bourdette. Interfacial transition zone in concrete. Advanced 
Cement Based Materials, 2:30-38, 1995. 

[23] I. Papayianni and M. Stefanidou. Strength-porosity relationships in lime-pozzolan mortars. 
Construction and Building Materials, 20:700-705, 2005. 

[24] B. Pichler and C. Hellmich. Upscaling quasi-brittle strength of cement paste and mortar: 
A multi-scale engineering mechanics model. Cement and Concrete Research, 41:467-476, 
2011. 

[25] B. Pichler, C. Hellmich, and J. Eberhardsteincr. Spherical and acicular representation of 
hydrates in a micromechanical model for cement paste: prediction of early-age elasticity 
and strength. Acta Mechanica, 203(3): 137-162, 2009. 

[26] K.L. Scrivener, A.K. Crumble, and P. Laugesen. The interfacial transition zone (ITZ) be- 
tween cement paste and aggregate in concrete. Interface Science, 12:411-421, 2004. 

[27] A. Sepulcre-Aguilar and F. Hernandez-Olivares. Assessment of phase formation in lime- 
based mortars with added metakaolin, Portland cement and sepiolite, for grouting of historic 
masonry. Cement and Concrete Research, 40:66-76, 2010. 

[28] M. Stefanidou and I. Papayianni. The role of aggregates on the structure and properties of 
lime mortars. Cement & Concrete Composites, 27:914-919, 2005. 

[29] J. Stransky, J. Vorel, J. Zeman, and M. Sejnoha. Mori-Tanaka based estimates of effective 
thermal conductivity of various engineering materials. Micromachines, 2:129-149, 2011, 
1101.4121. 

[30] E. Vejmelkova, M. Keppert, Z. Kersner, P. Rovnanfkova, and R. Cerny. Mechanical, 
fracture-mechanical, hydric, thermal, and durability properties of lime-metakaolin plas- 
ters for renovation of historical buildings. Construction and Building Materials, 31:22-28, 
2012. 

[31] A.L. Velosa, F. Rocha, and R. Veiga. Influence of chemical and mineralogical composition 
of metakaolin on mortar characteristics. Acta Geodynamica et Geomaterialia, 153:121-126, 
2009. 

[32] J. Vorel, V. Smilauer, and Z. Bittnar. Multiscale simulations of concrete mechanical tests. 
Journal of Computational and Applied Mathematics, 236:4882-4892, 2012. 

[33] V. Smilauer, P. Hlavacek, F. Skvara, R. Sulc, L. Kopecky, and J. Nemecek. Micromechan- 
ical multiscale model for alkali activation of fly ash and metakaolin. Journal of Materials 
Science, 46:6545-6555, 2011. 

[34] C.C. Yang. Effect of the transition zone on the elastic moduli of mortar. Cement and 
Concrete Research, 28:727-736, 1998. 

[35] A. Zaoui. Continuum micromechanics: Survey. Journal of Engineering Mechanics, 128:808- 
816, 2002. 



10 



A Herve-Zaoui solution 



The effect of coating on tlie meclianical properties enters tlie solution through the auxihary 
factors in Eq. (5), and and in Eq. (6). Here, these are provided in the closed form 
optimized for coding, utilizing the results and nomenclature by Herve and Zaoui [10]. Note 
that in order to keep the notation consistent, a^^^ corresponds to a property of the k-th phase, 
whereas a*^ denotes a quantity utilized in the Herve-Zaoui solution (independent of a^*^^). Also 
recall that we employ the local numbering of phases by the index i = [ii, 12,^3]^ introduced by 
Fig. 2. 

In particular, the volumetric part is expressed in terms of matrices 



Q2 = ^2q1 



with 



AT'' 



1 



3i^(»fe+i) + 4G'(»fc) 



for A; = 1,2. 



The matrices needed to evaluate the deviatoric part follow from 



p2 p2 
-'^11 -'^22 



p2 p2 ' 
-'^12-'^21 



^21 



p2 p2 
-'^11 -'^22 



p2 p2 ' 
-'^12-^^21 



where 



,fe-l 



p2 p2 
-'^22-f^ll 



p2 p2 ■ 
-'^12-'^21 



^21 



0] (A; = 1,2), 



p2 ^ m'^P\ 



The auxiliary matrix Ai"^ admits the expression: 



5(1 - i/fe+i)) 



3 




Mil 



with 



-20a'=(l -2i/(^'=+i)) 



i^fe)2(35fc _ 7cfc) 

5(1 - 2i/fe)) 
6'=(1 - 2z/(*fe+i)) 

5(1 - 2i/fe)) 
-i?(^'=)7(2a'= + UTa'^') 

70(1 - 2zy(*'=)) 
7Q,fc^fc)5(i - 2i/(*'=+i)) 

2(1 - 2i/fe)) 



M. 



23 7i?(*fc)7 

_ i?(*'')2(105(l - i/fe+i)) + 12a'=(7 - 10i/(^'=+i)) - 7e'=) 
34 — 



35(1 - 2zy(»fc)) 



-12a'= 

y 



4(/'' - 27a^') 
15i?fe)3(i_2z/fe)) 
-12a^(l - 2z/(^'=+i)) 
TR{ik)T{l - 2iy('k)) 

e'=(l - 2z^(*''-+i)) 
3(1 - 2i/(*fe)) 

-5a'' R^'"^^ {I - 2zy(*'=+i)) 



41 
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and 



a*^ = ^fa^^) (7 + 5Gfa))(7 - lOGfa+i)) - (7 - 10G(^*))(7 + 5G(*'=+i)), 

= y(7 + SG^^'^)) + 4(7 - 10G(^'=)), 
c*^ = (7 - 5G(^'=+i)) + 2(4 - SG^^^+i)) 



= (7 + 5G(*'=+i)) + 4(7 - 10G(''=+i)) 



6^= = 2(4 - 5G(^'=)) + y (7 - 5G(^'=)), 

Z*^ = (4 - 5G(^'=))(7 - 5G(^'=+i^) - y(4 - 5G(^'=+i))(7 - SG^^'^^), 
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